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1. INTRODUCTION 

Standard deviation and variance [1] measure how far a set of numbers are spread out from their 
average value and they give an indication of the fluctuation or chaotic nature of the numbers. For instance, 
let’s consider an array of numbers, Sı = [2, 2, 2, 2, 9, 9, 9, 9]. The standard deviation of the elements of S; is 
3.74 and the variance is 14. Now, let’s reorder the numbers of S; and let’s consider S2 = [2, 9, 2, 9, 2, 9, 2, 9]. 
For S2, the standard deviation is 3.74 and the variance is 14. Here, we see that standard deviation and 
variance do not change even if the numbers are reordered. The numbers in S2 has more changes with respect 
to its next number than in Sı, and Sz fluctuates more than Sı. However, standard deviation and variance 
cannot acknowledge this difference. 

Let’s consider a 2D case of numbers. An image is a 2D array of numbers i.e. pixel values [2]. Figure 
la and Figure 1b shows two artificial images where they have the same number of white and black pixels. In 
an 8 bit grayscale image, the white pixel has the value of 255 and the black pixel has the value of 0. Fig. 2a 
and Fig. 2b shows two natural grayscale images having different textures. The standard deviation, variance, 
and entropy of the image pixels in Figure 1 and in Figure 2 are shown in Table 1 and in Table 2 namely. 
From Table 1, we see that the standard deviation, variance, and entropy [3]-[5] of Figure la and Figure 1b are 
exactly the same. However, intuitively Figure 1b is more random or unpredictable than Figure la. From 
Table 2, we see that the standard deviation, variance, and entropy of Figure 2a and Figure 2b are very close 
even though they have different textures. For instance, Figure 3a and Figure 3b shows the changes of pixel 
values with respect to its adjacent left pixel value for row 127 for Figure 2a and Figure 2b namely. Here, we 
see that Figure 2b has more fluctuations of pixel values than Figure 2a. Quantitively, average absolute change 
for Figure 3a is 9.34 and Figure 3b is 22.49. 
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The standard statistical measures such as standard deviation, variance, and entropy cannot recognize 
the difference between the image pairs in Figure 1 and Figure 2 because they do not consider the order of the 
numbers. In an image, the order of the numbers (i.e. pixel values) are significant because changing their order 
will change the image. Thus if these standard statistical methods are used to find the randomness or 
fluctuation of the texture of an image, they may lead to misleading results. In this paper, an efficient solution 
of this problem is proposed. The proposed method sums up the changes of consecutive numbers and can 
recognize the different texture of images even though they have same standard deviation, variance, and 
entropy. The proposed method can be used to generate a feature vector in classifying textures of stone [6], 
wood [7], batik motif [8] etc. 


(a) 


Figure 1. Artificial images having the same number of white and black pixels. (a) Pixels distributed in an 
orderly fashion with fewer fluctuations; (b) Pixels distributed randomly having more fluctuations 


(b) 


Figure 2. Natural images having different textures. (a) Pixels distributed having fewer fluctuations; (b) Pixels 
distributed having more fluctuations 


Table 1. Standard Statistical Measurement of Figure 1 


Figure Standard Deviation Variance Entropy 
Figure la 127.50 16256.25 1 
Figure 1b 127.50 16256.25 1 


Table 2. Standard Statistical Measurement of Figure 2 


Figure Standard Deviation Variance Entropy 
Figure 2a 37.7375 1424.10 7.1794 
Figure 2b 37.7397 1424.30 7.2472 
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Figure 3. Changes of the pixel value at row 127 with respect to its left pixel. (a) Changes of pixel values for 
Figure 2a.; (b) Changes of pixel values for Figure 2b 


2. THE PROPOSED METHOD 

The proposed fluctuation measurement method, fin, for the 1D array is expressed in (1). Here x; is 
the im element of the array and N is the size of the array. The numerator contains the sum of the squared 
changes of consecutive numbers and the denominator contains the total number of changes. 


1 
(x, T Xin y 


ha EnA (1) 


For 2D case, such as an image, the proposed fluctuation measurement method, f2p, is expressed in 
(2), (3), and (4). Here (xj, yj) is the number at the in column and the jin row. X and Y are the total columns and 
total rows namely. Here in (2), fn represents the fluctuation in the horizontal direction and f, in (3) represents 
the fluctuation in the vertical direction. Finally, fop, in (4) is the mean of fn and f, multiplied by a factor n. For 
instance, n could be the maximum pixel value, i.e. 255 for an 8-bit image. 
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3. RESULTS 

Using (1), the fluctuation of S; is 7 and the fluctuation of S2 is 49 Thus, the proposed method can 
acknowledge the existence of the frequent changes in S2 and give higher fluctuation than S;, even though 
both S; and Sz have the same elements. Table 3 and Table 4 show the fluctuation for the image pairs in 
Figure 1 and Figure 2 using the proposed method namely. Table 3 shows that the fluctuation for Figure 1b is 
much larger than Figure la. From Table 4, we see that the fluctuations of Figure 2b is much larger than 
Figure 2a. Thus, the proposed method can recognize the different amount of fluctuations between the 
artificial image pairs in Figure 1 and between the natural image pairs in Figure 2. On the other hand - the 
standard deviation, variance, and entropy are unable to appreciate the different amount of fluctuations 
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between the image pairs in Figure 1 as shown in Table 1 and between the image pairs in Figure 2 as shown in 
Table 2. 


Table 3. Fluctuation measurement of Figure 1 using proposed method 


Figure Horizontal fluctuation (f) Vertical fluctuation (f) Fluctuation (fon) 
Figure la 0 1 127.50 
Figure 1b 4.87 5.80 1360.33 


Table 4. Fluctuation measurement of Figure 2 using proposed method 


Figure Horizontal fluctuation (fn) Vertical fluctuation (f) Fluctuation (fp) 
Figure 2a 41.05 33.39 9491.24 
Figure 2b 77.68 76.89 19707.70 


One existing method, mentioned in [9]-[11], for finding the fluctuation of an image could be to use a 
local standard deviation filter. This method returns a 2D array where each output pixel contains the standard 
deviation of the 3-by-3 neighborhood around the corresponding pixel in the input image. Symmetric padding 
is used for the pixels in the border. Taking an average of the returned 2D array can indicate the fluctuation of 
the image. As this method uses 3x3 neighboring pixels to calculate the standard deviation of each pixel, the 
pixel positions and orders are considered in this method. Local entropy filter can also implement in the same 
way where entropy is calculated instead of standard deviation. In Table 5, the fluctuation is measured in 
Figure 2a and Figure 2b using local standard deviation filtering, local entropy filtering and using the 
proposed method. Here, we see that all of these methods correctly reports higher fluctuation in Figure 2b than 
in Figure 2a. 

One shortcoming of the local standard deviation filtering and local entropy filtering methods are 
their time complexity. In these methods, the standard deviation (or entropy) of 9 pixels needs to be calculated 
for each pixel in the image. However, in the proposed method, only a difference between two pixels needs to 
be calculated. The local standard deviation filtering, local entropy filtering, and the proposed method were 
implemented (without implementing vectorization - for proper comparison) in a computer having Intel 
Pentium CPU 2117U running at 1.8GHz. The execution times are shown in Table 6. Here, we see that the 
proposed method is the fastest having the lowest execution time. The proposed method is approximately 8 
times faster than local standard deviation filtering and approximately 33 times faster than local entropy 
filtering method. Fourier spectra methods [2], [12] can be used to measure fluctuation, however, they have 
higher computational complexity than [9] and the proposed method - and will require more execution time. 


Table 5. Fluctuation measurement of Figure 2 


Figure Local Standard Deviation Filtering Local Entropy Filtering Proposed Method 
Figure 2a 10.71 2.82 9491.24 
Figure 2b 24.58 3.03 19707.70 


Table 6. Comparison of execution time (in Seconds) 


Figure Local Standard Deviation Filtering Local Entropy Filtering Proposed Method 
Figure 2a 3.39 13.50 0.39 
Figure 2b 3.38 13.56 0.37 


4. CONCLUSION 

Entropy, standard deviation, and variance may lead to misleading results if they are used to measure 
randomness or fluctuation. In this paper, a low complexity and fast fluctuation measurement method is 
proposed considering the order of the numbers for both 1D and 2D cases. The proposed method can be used 
in image processing applications such as classifying texture. 
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